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ABSTRACT 

Maps of the 450 /im and 850 /m dust continuum emission from three star-forming condensa- 
tions within the Lynds 1630 molecular cloud, made with the SCUBA bolometer array, reveal 
the presence of four new submillimetre sources, each of a few solar masses, two of which 
are probably Class I, two Class 0, as well as several sources whose existence was previously 
known. The sources are located in filaments and appear elongated when observed at 450 fim. 
They likely have dust temperatures in the range 10 to 20 K, in good agreement with previous 
ammonia temperature estimates. Attempts to fit their structures with power-law and Gaussian 
density distributions suggest that the central distribution is flatter than expected for a simple 
singular isothermal sphere. 

Although the statistics are poor, our results suggest that the ratio of 'protostellar core' 
mass to total virial mass may be similar for both large and small condensations. 



Key words: stars: formation - radio continuum: ISM - 
- ISM: individual: Orion B 



ISM: clouds -ISM: individual: L1630 



1 INTRODUCTION 

The Orion molecular clouds are the subject of extensive investi- 
gation on account of their proximity (400-500 pc) and the for- 
mation within them of stars of a variety of masses. An important 
unbiased survey of molecular condensations within the clouds is 
that of Lada, Bally and Stark (1991; hereafter LBS), who mapped 
the Orion B (L 1630) cloud with an angular resolution of 1.7 ar- 
cmin, observing emission from the J=2-l transition of interstellar 
CS. Excitation of this transition requires molecular hydrogen den- 
sities of order 10 4 — 10 5 cm -3 . Several of these condensations have 
themselves been mapped with substantially higher resolution and 
in molecular transitions requiring still higher densities for excita- 
tion (Zhou et al. 1991; Lis, Carlstrom & Phillips 1991; Chandler, 
Moore & Emerson 1992). 

Gibb & Heaton (1993; hereafter GH93) and Gibb et al. (1995; 
Paper I) mapped 6 of the LBS condensations in the J =3-2 or J=4- 
3 transitions of interstellar HCO + . These revealed an apparently fil- 
amentary structure containing embedded cores whose virial masses 
were typically a few solar masses. 

Paper I mapped the dust continuum emission from several of 
the cores in LBS23, while Gibb & Little (1998; hereafter GL98) 
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mapped the whole of LBS23 in J=2-l C O. Dust-derived masses 
for the cores were in good agreement with their virial masses de- 
rived from HCO + . The emission peaks evident in HCO + usually 
correlated well with those of the dust continuum but were faint or 
missing in C ls O. GL98 considered that these results implied deple- 
tion of the C 18 by factors of between 10-50. The lack of correla- 
tion between dust continuum and C 18 was similar to that found in 



the hotter source NGC2024 ( Mauersberger et al. 1992 ). The si mple 
interpretation as depletion in NGC2024 ( Mezger et al. 1992 ) has 
been questioned by Chandler & Carlstrom (19961 on the grounds 
that (a) the dust opacity law in the cores may be non-standard, 
(b) dust emission may become optically thick at the short submm 
wavelengths, but erroneously interpreted as thin, (c) the cores may 
contain unresolved components which are highly optically thick in 
C ls O but interpreted as thin so the true molecular abundance is 
much higher. 

The recent availability of the Submillimetre Common User 
Bolometer Array (SCUBA: Holland et al. 1999) permits rapid map- 
ping of dust emission from the LBS condensations with greatly im- 
proved sensitivity. This means that large areas within them can be 
mapped to determine the gross dust distribution and further proto- 
stellar cores may be readily identified. We here describe the results 
of a pilot study in which selected areas of our LBS condensations 
have been mapped with SCUBA. 

The motivation for this study was: 

(a) to search for dust components near, for example, outflow 
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Figure 1. SCUBA maps of LBS23S at 450 fim (left) and 850 fim (right). The known sources HH25MMS, SSV59 and HH26IR are marked by crosses, while 
the components of HH25 and HH26 from Davis et al. (1997) are marked by open squares. The greyscale is in Jy/beam and is shown as a bar across the top of 
each figure, extending from — 3<r to the maximum (4.0 Jy/beam at 450 fim and 1.4 Jy/beam at 850 fim). Contours are -3, 3, 5, 7, 9, 12, 15, 18, 21, 24 (450 
lim) and 5, 10, 15, 20, 25, 30, 35, 40 (850 fim) times the noise level of 150 mjy per beam and 36 mjy per beam respectively. 



sources, near-infrared sources, and Herbig-Haro objects, where 
they have not previously been detected, and for new protostellar 
objects. 

(b) to compare the observed dust emission with mod els of the 
structure of pr otoste llar cores such as those of Shu (1977), or 
Crutcher et al. ( 1994 ). The simultaneous mapping of the cores at 
850 and 450 fims facilitates some separation of dust mass and tem- 
perature. 

(c) to see whether the morphology of large scale filaments frag- 
menting into protostellar cores evident in e.g. NGC2024 and 
LBS23 extends to smaller condensations. 

(d) to compare the numbers of protostellar cores in large and 
small condensations in the same molecular cloud and to see if stars 
form preferentially in the larger condensations. 

(e) to map fields for which detailed HCO + and C ls O maps exist 
already, to determine molecular depletion in different source com- 
ponents. A detailed analysis will follow in a later paper. 



2 OBSERVATIONS 

The observations were made using the submillimetre bolometer ar- 



ray receiver SCUBA (Holland et al. 1999) operated in service mode 



at the James Clerk Maxwell Telescope over the period 24 to 26 Au- 
gust 1997 and 20 and 23 September 1997. The 32 element 850 /im 
and 96 element 450 fim arrays were used in 'jiggle-mode' to make 
fully-sampled maps of 2.2-arcmin diameter-fields in 6 of the LBS 
cores. The chop throw was 120 arcsec in azimuth. The flux scales 
were established by observations of Uranus at the start of each ob- 



serving session and skydips were used to calculate the zenith opti- 
cal depth at the start and end of each session. 

The night of 26 August was very humid and calibration pro- 
duced results that differed by of order 50 per cent (lower) com- 
pared with the previous and following nights. The weather on the 
nights of 20 and 23 September was poor and the 450 fim data could 
not be used. Thus the images presented were constructed from data 
recorded on the nights of 24 and 25 August 1997 only. On these two 
nights, the estimated calibration uncertainty is of order 10 per cent 
at 850 fim and 30 per cent at 450 fim. The beam of the telescope 
was determined to have full-width at half-maximum (FWHM) di- 
ameter of 15.0±0.5 arcsec at 850 fim and 8.1 ±0.2 arcsec at 450 
^m. The total time on each source was of order one hour, except 
for LBS7, 13 and 16. Due to poor weather the time on source for 
these fields was only 21 minutes. 

The typical noise level at 850 fim and 450 ^im was 34 mjy per 
beam and 150 mjy per beam respectively for LBS17, 18 and the 
two fields in LBS23. This translates to a 3<r mass sensitivity of 0.6 
M0 for temperatures of 10 K, assu ming that the dust emissivity is 
as described by Hildebrand (1983). The noise level for LBS7, 13 
and 16 was higher at ~80 mjy per beam at 850 fim; no emission 
was detected in these fields. 

Table [l] summarises the fluxes and positions of the dust emis- 
sion sources detected in the LBS cores. 
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Figure 2. SCUBA maps of LBS23N at 450 /Um (left) and 850 fim (right). The known sources HH24MMS, SSV63E and SSV63W are marked by crosses while 
the Herbig-Haro objects HH24A-D are marked by open squares. The greyscale is in Jy/beam and is shown as a bar across the top of each figure, extending 
from -3<r to the maximum (5.3 Jy/beam at 450 fan and 2.5 Jy/beam at 850 fan). Contours are -3, 3, 5, 7, 9, 1 1, 15, 20, 25, 30 (450 /im) and 10, 15, 20, 25, 
30, 40, 50, 60, 80, 100 (850 fj,m) times the noise level of 160 mjy per beam and 30 mJy per beam respectively. 
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Figure 3. SCUBA maps of LBS 17 at 450 fim (left) and 850 fim (right). The water maser of Haschick et al. (1983) is marked by a cross. The greyscale is 
in Jy/beam and is shown as a bar across the top of each figure, extending from — 3cr to the maximum (3.0 Jy/beam at 450 fim and 1.4 Jy/beam at 850 fim). 
Contours are -3, 3, 5, 7, 9, 12, 15, 18, 21, 24 (450 /im) and -3, 5, 10, 15, 20, 25, 30, 35, 40 (850 (im) times the noise level of 130 mjy per beam and 35 mjy 
per beam respectively. 
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Figure 4. SCUBA maps of LBS18 at 450 fim (left) and 850 fim (right). The 3.5-cm source detected by Gibb (1999) is marked by a cross. The grayscale is 
in Jy/beam and is shown as a bar across the top of each figure, extending from — 3cr to the maximum (1.7 Jy/beam at 450 fim and 0.91 Jy/beam at 850 fim). 
Contours are -3, 3, 5, 7, 9, 12, 15, 18, 21 (450 ^m) and 5, 8, 1 1, 14, 17, 20, 23, 26 (850 fim) times the noise level of 150 mjy per beam and 33 mjy per beam 
respectively. 



3 RESULTS 

Figs. 1 to 4 show the SCUBA maps at 450 and 850 for the 
four LBS fields in which emission was detected. Table |l] lists the 
parameters derived from two-dimensional gaussian fits (using the 
task IMFIT in AIPS) for the sources detected in the LBS 17, 18 and 
23 fields. The source dimensions have not been deconvolved. The 
(formal) uncertainties in the positions and dimensions arising from 
the fitting procedure are of order 0.5 arcsec. The main features of 
each source will now be discussed in turn. 



3.1 LBS23S 

This field i s centred on the class source HH25MMS (G ibb & 
Davis 19 98) and contains th e class I infrared sources SSV59 and 
HH26IR ( [Davis et al. 199% . The purpose of mapping this region 
was to compare the SCUBA maps of HH25MMS (Fig. 1) with 



the existing UKT14 map (Paper I) and also to search for the driv- 
ing source of the southern bipolar outflow in HH24-26 (GH93). 
The morphology of HH25MMS is consistent in both maps, and the 
maximum l ies w ithin a couple of arcsec of the 3.4-mm peak of 
Choi et al. (1999). Two further weaker components are evident in 
Fig. 1; they are associated with HH26IR, which is the source of the 
southern outfl ow (GH93), and SSV 59, which also drives a molec- 
ular outflow (Gibb & Davis 199"). Edge effects make it difficult 



to make accurate measurements of the flux density of HH26IR as 
we have not observed the entire source. Recent 1300-^m and 350- 
fim observations of the whole of LBS23 by Lis, Menten & Zylka 
(1999) reveal a structure consistent with what we observe in Figs. 1 
and 2. They show that HH26IR also appears to be embedded within 
a clump elongated perpendicular to the CO flow. 

A curious point is that the radio continuum source and origin 
for the HH25MMS molecular outflow is not coincident with the 
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submillimetre continuum peak, being offset by ~6 arcsec in dec- 
lination, somewhat larger than the mean pointing error (3 arcsec). 
If genuine, this observation may indicate the presence of a second 
submillimetre source within LBS23S, which may be younger than 
HH25MMS (on account of there being no observable outflow). Al- 
ternatively, it may be a hot-spot heated by the impact of the out- 
flow from HH25MMS. Clearly higher resolution observations are 
needed to clarify this issue. Lis et al. (1999) observe this in their 
1300-/im map, but it seems less pronounced in their 350- /im data. 



3.2 LBS23N 

This field is centred on the HH24 complex, and includes the class 
protostar HH24MMS. The two components HH2 4MMS and 
SSV63W detected at 1300 fxm by Chini et al. (Il993|) are clearly 



mapped (Fig. 2). There is good agreement between the 850 /j,m 
SCUBA and 800 /im UKT14 (Paper I) peak fluxes for the unre- 
solved source HH24MMS. The 450 /im map reveals the source to 
be extended in a N-S direction. While the flux densities are also 
in good a greement with the UKT14 data of Ward-Thompson et 
al. (1995) the SCUBA 450 /im observations show the source as 
much more compact. We note that the direction of th e com pact jet 
discovered by Bontemps, Ward-Thompson & Andre (1990 is nei- 
ther parallel nor perpendicular to the elongation at 450 /im. A new 
weak source is seen 40 arcsec north-west of SSV63 which shows 
up weakly at 450 /im as well. We label this clump HH24NW (Ta- 
ble [l]) although its parameters are not particularly well determined. 
Lis et al. (1999) also detect this source. 



3.3 LBS17 

Fig. 3 presents the submillimetre map of the south-eastern portion 
of the HCO + ridge of Paper I. The J=3-2 HCO + map produced a 
peak largely missing in C 18 (Paper I; Gibb & Little 2000). How- 
ever, this peak stands out as a dust source in the SCUBA map. 
There is a bipolar outflow associated with this source, observable 
in HCO+ and CO (Gibb & Little 2000). The elongation of the dust 
source, observable at 450 /im, is perpendicular to the direction of 
the bipolar outflow. No IRAS, optical or near-infrared emission ap- 
pears to be associated with the source - so it is likely to be young, 
possibly protostellar. Gibb & Little (2000) showed that it satisfies 
the criteria for being a class candidate. We refer to this source 
as LBS17H, to distinguish it from another source detected in the 
LBS 17 condensatio n to th e north of the area covered by our map 
by Launhardt et al. ( 199( ). This second source lies near a peak of 
our HCO + map (Paper I). A ridge of dust emission (Fig. 3) is per- 



ceived to run NW towards an extended plateau leading on to this 
peak, a feature also evident in HCO + . A poorly-defined peak is 
seen in the northwest of the 850- /im image although this may be an 
artefact of the mapping procedure. However, it is coincident with 
one of the HCO + clumps of Paper I, clump B6 although in keeping 
with the revised nomenclature we label it LBS17F in Table |l|. A 
more extensive map is required to determine its true parameters. 



3.4 LBS18 

LBS18 has one of the smallest virial masses of the LBS cores (14 
M©). In dust emission (Fig. 4), as in J=3-2 HCO + (Paper I) and 
J=2-l C 18 (Little et al., in prep), the source has a 'cometary' 
appearance, showing a N-S elongation with a co mpact component 
at one end, which could be protostellar. Gibb ( 1999| ) detected a 



Table 2. Masses derived from simple isothermal analysis assuming a tem- 
perature estimated from previous NH3 data and an emissivity index of 2. 
The mass within 8 arcsec was calculated from the 450 /im flux density; 
that within 15 arc sec fro m the 850 /tm flux density. The Class label is that 
define d by Lada (1991) and updated by Andre, Ward-Thompson & Bar- 
sony (19931. 



Source 


'id 
(K) 


M(M ) 
(8 arcsec) 


M(M ) 
(15 arcsec) 


Class 


Outflow 


LBS17H 


15 


0.9 


2.1 





Y 


LBS18S 


15 


0.5 


1.4 


9 


N? 


HH25MMS 


18 


0.8 


1.6 





Y 


SSV59 


18 


0.2 


0.5 


I 


Y 


HH24MMS 


12 


2.9 


7.9 





Y 


SSV63 


12 


1.3 


2.3 


I 


Y 



radio continuum source within this core using the VLA but it lies 
within the main north-south ridge, and is not coincident with the 
compact component. 



4 ANALYSIS 

4.1 Simplest model - spherically-symmetric, isothermal 
cores 

To determine the mass and temperature of the clumps requires a 
knowledge of the emissivity law, and that the distribution of emis- 
sion should be well resolved and observed with high signal-to-noise 
ratio at at least two different wavelengths. These conditions are not 
satisfied for our observations. The emissivity law is uncertain and 
even if this were not the case uncertainty in calibration (dominated 
by the 450 /im measurements) would be a strongly limiting factor 
in fixing the temperature and hence the mass. The clumps are barely 
resolved and the temperature within them may well vary with posi- 
tion, so that only an oddly weighted beam average can be derived. 
How then are we to proceed? 

Initially we have assumed the clumps possess the simplest 
possible structure. They were assumed to be optically thin and at 
a constant temperature, which we obtain from previous ammonia 
observations (Menten, Walmsley & Mauersberger 1987; Verdes- 
Montenegro & Ho 1996). Where no prior temperature data exist 
we have assumed an intermediate value of 15 K. Better estimates 
of the dust temperature may be obtained from g reybod y fitting if 
enough data exist. We employed the Hildebrand (1983) emissivity 
value of 0. 1 cm 2 g~ 1 at 250 /im to estimate the appropriate values at 
850 /im and 450 /im assuming a value of (5 of 1.5, the median value 
typical of dust emission in dense cloud cores (e.g. Huard, Sandell & 
Weintraub 1999; Chandler & Richer 2000). The core mass is given 
by F v d 2 /{kvB v ^t) where F„ is the flux density at frequency v, d 
is the distance to the source, k„ is the absorption coefficient and 
B V .T is the Planck function for a blackbody of temperature T. 

Even if these assumptions closely model reality, the calibra- 
tion error alone limits the determination of mass to within 50 to 
100 per cent. The assumed temperatures and derived masses are 
shown in Table ^| Despite the uncertainties, it seems that the dust 
cores probably have a mass within a 0.03 pc diameter region of or- 
der 1-2 Mq . The mass derived for HH24MMS seems anomalously 
high, which led Chini et al. (1993) to suggest that the absorption 
coefficient is larger towards this source (possibly due to a lower 
value for /3 or a coagulation of dust grains). 

We have made greybody fits to the spectral energy distribu- 
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Table 1. Source parameters, including identification where known. The source dimensions were derived from 2-dimensional gaussian fits and list the FWHM 
of the major and minor axes (in arcsec) and the position angle of the major axis in degrees east of north. Angles are rounded to the nearest integer. The last 
three columns list the flux density (in Jy) within an aperture of either 8 or 15 arcsec. 
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PA 
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) 


core 




RA(1950.0), Dec(1950.0) 


arcsec 2 


degrees 


arcsec 2 


degrees 


450 (im 


850 /im 


450 /J. 


17 


LBS17H 


05 h 43 m 57.1 s , -00°03'44.3" 


18x16 


26±2 


12x9 


49±3 


3.0 


1.4 


5.4 


17 


LBS17F 


05 h 43 m 54.2 s , -00°02'48.2" 


32x23 


-27±3 








0.4 




18 


LBS18S 


05 h 43 m 54.2 s , +00° 18'23.1" 


26x20 


22±4 


14x10 


31 ±3 


1.7 


0.9 


3.5 


23 


HH25MMS 


05 h 43 m 33.8 s , -00° 14'46.5" 


29x15 


-8±1 


15x13 


-1±9 


3.7 


1.4 


S.7 


23 


SSV59 


05 h 43 m 31.4 s ,-00°15'27.3" 


25x22 


13±6 


11x10 


79±43 


1.1 


0.4 


2.2 


23 


HH26IR 


05 h 43 m 30.1 s ,-00°16'00.l" 


27x17 


42±2 


9x8 


-57±64 


1.4 


0.6 




23 


HH24MMS 


05 h 43 m 34.7 s , -00° 1 1'48.8" 


17x16 


2±5 


13x9 


9±2 


5.1 


3.5 


9.6 


23 


SSV63 


05 h 43 m 34.4 s , -00° 1 1'07.3" 


23x17 


63±5 


16x10 


79±3 


2.4 


1.0 


5.1 


23 


HH24NW 


05 h 43 m 32.3 s , -00° 10' 43.1" 


24x21 


21±7 








0.5 





Table 3. Best fit parameters for Gaussian (columns 2—7) and power-law (columns 8-11) oblate models for each object modelled. Note that the mass estimates 
have been corrected for helium by multiplying the H2 results by 1.25. For the Gaussian model c n represents the maximum density and d n and e n the radii (to 
1/e) on the major and minor axes. For comparison, the HWHM beam radius at 450 micron is 0.008pc.The i z and i T columns show the indices for the power 
law njj 2 distribution in each direction (i.e. when 9 = w and = respectively in equation ^ ). c n is then the density at the inner cut-off radius O.OOlpc. The 
temperature distribution has been assumed to be spherically symmetric and varies according to Tk = axR? 3 + bxR~ 03 where R is the distance to the 
centre of the cloud measured in parsecs, ay being zero for the power law models. 



Object 


Mass 


Density 


Scale 


Temperati 
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Density 


Index 




Temperatu 




(M ) 


c n (cm -3 ) 


d n (pc) 


e n (pc) 


aj- 


br 


c n (cm -3 ) 




i r 


bx 


LBS17H 


16 


2.4xl0 7 


0.014 


0.010 


10 


1.6 


1.1x10 s 


-1.4 - 


■1.2 


2.2 


LBS18 


21 


9.8xl0 6 


0.025 


0.010 





2.1 


2.3 xlO 8 


-1.7 - 


■1.3 


1.7 


HH24MMS 


50 


2.6xl0 8 


0.009 


0.007 


15 


1.0 


5.6x10 s 


-1.6 - 


■1.2 


1.9 


HH25MMS 


12 


1.4xl0 6 


0.045 


0.018 


5 


4.3 


1.3xl0 7 


-1.1 - 


-0.7 


3.9 


SSV63 


8 


1.9xl0 6 


0.030 


0.015 


10 


2.9 


3.0xl0 7 


-1.6 - 


■1.3 


3.0 



tions using our SCUBA and IRAS HiRes-processed data to place an 
upper limit on the dust temperature (see Gibb & Little 2000, Gibb 
& Davis 1998 for greybody fits to two of our SCUBA sources). In 
general this temperature lies between 20 and 30 K, depending on 
the value of f3 for the fit, with the clumps containing IRAS point 
sources tending to have higher upper limits than those dust clumps 
with no IRAS point source. Unfortunately the greybody fits can- 
not provide more accurate estimates of the temperature because the 
IRAS fluxes are measured with a beam considerably larger than 
that of our JCMT observations, one which often encompasses more 
than one of our SCUBA clumps. 

If a power law variation of density with radius were assumed 
and the temperature constant this would imply a power-law index 
of -1.7, a value intermediate between the value expected for a sin- 
gular isothermal sphere in hydrostatic equilibrium (-2) and one in 
free-fall collapse (-1.5). If temperature increases inward the power 
law would be flatter. If a Gaussian distribution is assumed, an an- 
gular diameter for the sources ~12 arcsec (0.023 pc) is indicated, 
probably too close to the SCUBA beamsize to be meaningful. 



4.2 More detailed modelling of the emission 

Although the data indicate that the observations have only partially 
resolved most of the sources (particularly in the narrow direction) 
the existence of simultaneously observed well registered maps at 
both 850 and 450 fans justified more detailed modelling which 
would be capable of taking into account density and temperature 
gradients and finite optical depth. In particular we wished to com- 
pare the predictions of both Gaussian and power-law models, as 



they represent popular paradigms for the structure of such objects. 
We might attempt to distinguish between power law and Gaussian 
distributions by fitting the 850 fan and 450 fim observations simul- 
taneously. 

Accordingly, we have written a program that computes fluxes 
for model clouds with axial symmetry, described in the Appendix. 
To compare the model output with the telescope data two cuts were 
made on each of the SCUBA maps, one down the long axis of the 
cloud and one perpendicular to this. These cuts were then repeated 
on the model output and the results compared. The best fit oblate 
Gaussian and power law models for each object modelled are pre- 
sented in Table ^. Given the uncertainties in the calibration, which 
will directly propagate into the magnitude of the density and tem- 
perature in the model, we emphasize that in this section we seek 
to explore the possible forms of the density variation, rather than 
absolute values. 

We tried both the Hildebrand (13=2) and Testi and Sargent 
(1998) (/3=1.1) versions of the emissivity function (see Appendix) 
to see whether the choice of j3 affects our results significantly. For 
the Gaussian distribution there is no significant observable change 
in the any of the output cross-sections. However, it is necessary to 
decrease the density by a factor of approximately 2 and increase 
temperature by a factor of approximately 1.5 in order to retain the 
correct intensities which has the result of decreasing the mass of 
the cloud by around 50 per cent. Therefore there does not seem to 
be any clear way of distinguishing between different values of /3 
using this data set. 

We have also investigated the differences caused by assuming 
an oblate or a prolate spheroid as the cloud shape. Unfortunately 
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there does not seem to be any way of distinguishing between the 
two models as the best fit line profiles are virtually identical in both 
cases. However, the parameters needed to obtain the best fit are sig- 
nificantly different. The Gaussian prolate case needs a higher den- 
sity by a factor of between 1.3 and 3 (the more elliptical the cloud, 
the larger the difference). This in turn causes a higher cloud mass 
(up to 15 per cent). There is also a corresponding slight decrease in 
the temperatures required. All other parameters remain the same. 
For the power law distributions the same comments apply except 
that the increase in density caused a larger increase in the cloud 
mass (up to 50 per cent). 



4.2.1 Features of the modelling 

The results of the model fits are given in Table ^. In attempting the 
fits it was found that: 

(1) At least two of the sources (LBS18S and HH25MMS) 
showed evidence for two or more components along their major 
axes. 

(2) There was usually evidence of extended 'ridge' or 'wing' 
emission (up to 25 per cent of the peak) which could not be ac- 
counted for in the modelling carried out here. 

(3) Gaussian fits are usually slightly better than power laws in 
one significant respect, described below. 

Fig. ^ shows the fits for LBS 17. The presence of a ridge con- 
necting the main clump to a secondary clump to the north west 
causes an extended low level emission wing in the 850 ^tm i-axis 
cut which no effort has been made to fit. Ignoring this, the central 
part of the a>axis cut and the whole of the t/-axis cut is almost per- 
fectly fit by a Gaussian distribution profile. By contrast we were 
unable to find a power law fit that successfully modelled both the 
450 fim and 850 fim maps simultaneously. 

For power laws the 850 /im y-axis cut cannot be made narrow 
enough without causing the 450 fim ?/-axis cut to become too nar- 
row. The smaller 450 beam weights the central material more 
heavily - too heavily for the inwardly increasing power law. 

Of the five sources which were modelled this feature was ob- 
served along the major axes of SSV63, LBS17H, and LBS18S, 
and the minor axes of HH25MMS, HH24MMS, and LBS17H. Al- 
though the effect is not large it is noticeable and we believe it to 
be a significant discriminant. This would suggest that the density 
law flattens towards the centre of the clump (as in a Gaussian), 
although thus far such a requirement has only been determined 
for pre-collapse objects (e.g. Ward-Thompson et al. 1994). Future 
higher-resolution observations are clearly necessary to determine 
the density variation in the inner regions of these clumps. 

Investigation on dummy data showed that the effect of fitting 
a single component model to a two component source (i.e. main 
component plus 'ridge') was to produce a behaviour in the opposite 
sense to that observed. The one object where a power law seems to 
fit reasonably well (SSV63) has a less pronounced peak on the map 
(Fig. ^J) than the others. 



5 DISCUSSION 

5.1 Dust temperatures 

From simple isothermal greybody modelling the derived dust 
temperatures in LBS23 are generally consistent with those de- 
rived from HCO + and NH3 observations except in the case of 
HH24MMS, which yields a low dust temperature (Table ^). As 



mentioned above, Chini et al. (1983) attribute this to an enhanced 
absorption coefficient. However, it should be noted that greybody 
fitting of the spectral energy distribution yields more accurate re- 
sults as with enough data, the emissivity and its frequency depen- 
dence may be constrained in addition to the temperature (e.g. Ward- 
Thompson et al. 1995). 

The dust temperatures tend to be rather lower than the CO- 
derived temperatures, which, on account of the very high CO op- 
tical depth, are likely to refer to the outer parts of the cloud. Al- 
though external heating of these outer parts therefore seems likely, 
it should be noted that the temperatures estimated for the ridges of 
LBS17 and LBS18 show little or no evidence of it. This fact in- 
dicates that the majority of the dust emission arises from the cold 
interior of the cores observed here. 



5.2 Comparison with protostellar models 



The canonical protostellar collapse model of Shu (1977) predicts 
inside-out collapse in an isothermal spherical envelope, leading to 
a density distribution with a power law of index -1.5 in the col- 
lapse region while outside it the index has the isothermal value of 
-2.0. The effect of magnetic fields and/or rotation will obviously 
complicate this simp le picture. For example, the simulations of 
Crutcher et al. (1994), which include the effects of magnetic fields 



and ambipolar diffusion, indicate modifications of these laws to - 
1.3 and -1.7. Flattening would be expected leading to elongation in 
projection onto the plane of the sky. 

The observations of the clumps do indeed reveal elonga- 
tion which is particularly evident at 450 /im where the 8 arcsec 
beamwidth of the telescope is equivalent to 3200 AU. However, 
HH25MMS, HH24MMS and LBS 18 are elongated in the direction 
of the filamentary structure in which they are embedded, suggest- 
ing perhaps that they might be prolate. In SSV63 and LBS17H the 
elongation appears to be perpendicular to the filaments. 

When the observable extent of a source is not much greater 
than the beamwidth it is perhaps invidious to seek to distinguish be- 
tween Gaussians and power laws. The central regions of the clumps 
can be fitted by Gaussians. It is hard to separate wing emission from 
a clump from emission due to the embedding filament. It is clear 
that, if the densities are interpreted as power law distributions, the 
power law indices vary with direction. The indices lie in the range 
-0.7 to -1.2 in the direction of elongation and -1.1 to -1.7 perpen- 
dicular to it. These values have been derived assuming an r -0 ' 3 
temperature dependence. If the clumps were taken to be isother- 
mal their magnitude would increase by about 0.3. Given that in the 
central regions Gaussians are somewhat favoured anyway, the least 
that can be said is that the results do not agree very well with the 
simple Shu model. 

Questions pertinent to our results include: (a) are the cores 
oblate or prolate?, and (b) what has led to their orientation relative 
to the filaments in which they are embedded? While a number of 
papers have been published examining the fragmentation of fila- 
ments (e.g. Fiege & Pudritz 2000; Nakajima & Hanawa 1996), few 
have attempted to analyze the structure of the fragments and in par- 
ticular look at the issue of prolate or oblate cores. However, Bon- 
nell, Bate & Price (1996; hereafter BBP) have performed computer 
simulations of the dynamic evolution of initially prolate, magnetic 
field- and rotation-free filaments into protostars, which may be use- 
fully related to our observations. The evolution depends primarily 
on the ratios of the Jeans radius to the minor and major axes of the 
filament, and they show that under certain conditions it is possible 
to produce an oblate core embedded within a prolate condensation. 



© 2000 RAS, MNRAS 000, \l 



8 Robin R. Phillips, Andy G. Gibb and Leslie T. Little 



ntensity along X — axis of the 450 cloud 



ntensity along Y— axis of the 450 cloud 




-10 10 

Position (") 



20 





I I I I 

\ 


— 
- 


- : 

r 


\ 


- 
- 


1 : 


■ V 




7 ' 






- // ■' 


■ V 




/ // 






/ // 






/ /•' 


A 


- 


/ / 




\ - 


/ 




\ N 










-20 -10 10 

Position (") 



20 



ntensity along X-axis of the 850 cloud 



ntensity along Y— axis of the 850 cloud 




-10 10 

Position (") 



20 



I I I I I I I 














\ 












\ 










1 


\V. 




/// 


V 




■'// 






• 7 






••'// 






■ 1 






/I 




\ ■■ 


' / / 




V 


- .•' / / 




V ■ - 


/ / 




\\ 


/ / 




\\ - 


? / 





II 



-20 -10 10 

Position (") 



20 



Figure 5. Cuts along the major and minor axes of LBS17H. The solid and dotted lines are respectively the modelled results from a Gaussian density distribution 
and a power law distribution. The dashed lines are cuts taken from Fig. j?| The strips are the la and 3cr noise levels (from Fig. ^). The beamsizes used were 8 
arcsec for 450 /Um and 15 arcsec for 850 fim. 



BBP give arguments to suggest that rotation would be unlikely 
to modify hugely these results. However, despite the claims of BBP, 
it is likely that magnetic fields may do so although their effects 
would depend on the field strength, direction and the Alfven speed. 
For example, if the field lines were perpendicular to the axis of the 
filament then collapse across the field lines towards a single object 
might be inhibited on dynamical timescales, resulting in a string of 
individual cores which each form a protostar. 

Can these theoretical considerations be related to our observa- 
tions? First, we note that the effect of projecting both prolate and 
oblate object onto the plane of the sky is to decrease their apparent 
asymmetry. Both the filaments and the envelopes are, on average, 
more asymmetric than they appear. Second, the typical core size 
of a few x 10 3 au is similar to the scale of the flattened structures 
in BBP's simulations. Thirdly, like LBS23 and NGC2024, LBS 17 
and LBS 18 both have an elongated filamentary structure with pro- 
tostellar cores appearing to form as a result of fragmentation along 
their length. Within the framework of BBP's simulations, their orig- 
inal filamentary structure was sufficiently massive that fragmenta- 



tion has occurred along their length, as may have also happened in 
LBS 18 (see Fig. 4). One would therefore expect that the condensa- 
tions in these objects should be regarded as BBP's 'intermediate' 
prolate objects aligned with the original filaments. 

HH25MMS, HH24MMS and LBS18S are elongated in the di- 
rection of the filamentary structure in which they are embedded, 
suggesting perhaps that they might indeed be prolate 'intermediate' 
cores which form after fragmentation of a prolate filament (BBP). 
But in SSV63 and LBS 17H the elongation appears to be perpendic- 
ular to the filaments - within the framework of BBP's analysis they 
would be oblate 'ultimate' cores. Both these cores are observed 
at the ends, rather than centres, of filaments, however. Could it be 
that they have 'swallowed' all the material on one side of them? 
However, LBS18S is also observed at the end of a filament, but is 
elongated along the filament. 

As there is evidence of star formation all along LBS23 it 
seems likely that protostars have sufficient time to form within 
the intermediate envelopes. Their mutual attraction could lead fi- 
nally to a stellar cluster - not BBP's merging into a single object 
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which results from the symmetry of their initial starting conditions. 
BBP also note that the initial conditions require the accumulation 
of many Jeans masses in a time much shorter than the dynamic 
time scale, which might be achieved, for example, by cloud-cloud 
collisions, a s prop osed for the formation of the filament TMC 1 by 
Little et al. ((1978 



6 CONCLUSIONS : SEARCHING FOR PROTOSTARS 

Searching LBS cores with SCUBA is a powerful method for detect- 
ing protostellar candidates in the Orion clouds. Several new ones 
have been discovered in very little observing time. On the sim- 
plest interpretation the clumps observed have masses (in the 15 
arcsec beam) varying from a few (SSV59, SSV63, HH26) to 50 
(HH24MMS) solar masses. Of the newly detected sources SSV59 
and HH26IR are probably Class I objects, while LBS17H is a new 
Class source. LBS 1 8S has no near-infrared, IRAS or radio source, 
nor does it appear to have an outflow and so it is very young, per- 
haps pre-protostellar. Further studies of these objects, e.g. by inter- 
ferometry, are required to confirm their evolutionary status. 

One (LBS 18) out of the 4 small LBS cores shows evidence 
of a protostellar clump of several solar masses. There are two in 
LBS17 (including the 1300 fim source of Launhardt et al. 1996) 
and several in LBS23. On the basis of present limited observations 
(4 small LBS cores, 2 large LBS cores) the ratio of 'protostellar' 
clump mass to the total virial mass of the filament does seem not 
very different for large or small cores - but the statistics are not 
good. A proper census of the LBS cores, in which both small and 
large cores are mapped in their entirety, feasible with SCUBA, will 
allow a mass spectrum of these clumps to be determined. The next 
obvious step is to carry out a wide field, unbiased survey of L 1630 
to determine the mass distribution on 15-arcsec scales. 
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APPENDIX A: DUST RADIATIVE TRANSFER MODEL 

Our radiative transfer model is a simple modification to the ax- 
isymmetric molecular line transfer code of Phillips & Little (2000) 
to cope with dust continuum emission. 

The emissivity, E v , of a small element i is written as 



2hu s 



c 2 (e J 



(Al) 



where the dust is at a temperature The cloud to be modelled 
is split into geometrical subsections. The geometrical subsections 
chosen for the cloud are a series of stacked concentric rings. Each 
subsection is assigned a temperature and a molecular hydrogen 
density according to the desired model. For any desired line of 
sight the positions of the intersections of the line of sight with the 
boundaries of the subsections are then calculated. It is then possible 
to integrate the emission along the line of sight simply by adding 
up the contributions from each segment on the line of sight (Ah) 
according to 



E u ,iAke 



-V* ■ 



(A2) 



© 2000 RAS, MNRAS 000, [l|-| 



10 Robin R. Phillips, Andy G. Gibb and Leslie T. Little 



where E Ul i and T u ,j are the emissivities and optical depths for seg- 
ment i (or j) respectively. 

The beam-averaged emission is computed by taking a grid of 
positions over the beam, calculating the emission from each and 
summing the results with appropriate weightings. For all the mod- 
elling that follows the beam is assumed to be a perfect Gaussian. 
Whilst this is a good assumption for the 850 ftm data the JCMT 
beam deviates from a perfect Gaussian at 450 due to imper- 
fections in the dish. Observations of Uranus showed that the beam 
had a substantial 'error' component containing ~50 per cent of the 
power and spread over a 50 arcsec-diameter region. However, al- 
though the precise shapes of the 450 fim cuts may be slightly inac- 
curate (mostly at the lower levels), the difference is almost certainly 
not enough to significantly affect our results. 

Two different models for the molecular hydrogen distribution 
were tested. The first is a power law model of the form 



"Ilo 



c n D~ 



where D=R/Ri n , 6 = tan 



)1 (A3) 

j and r and z are the cylindrical 
co-ordinates of the position of interest in the cloud with the z-axis 
being the axis of symmetry in the cloud and r the distance from the 
axis, while R refers to the distance from the cloud centre measured 
in parsecs. The z axis is taken to lie in the plane of the sky. Choos- 
ing appropriate values of d n and e n enables the power law depen- 
dence to vary with direction, thus allowing an elongated cloud to be 
described. By changing the sign of e it is possible to model either 
a prolate or an oblate cloud. Although such a description has the 
major advantage of enabling easy comparisons with various pre- 
dictions for the power law in each direction it is difficult to obtain 
smooth elliptical contours of equal density. It is necessary to give 
the value of c„ for a position, R[ n , close to the centre of the cloud 
rather than at the outer edge as many authors have done (e.g. Lit- 
tle et al. 1994) since specifying the density at the outer edge causes 
the two different power laws to give two different density values 
at the cloud centre. The apparently odd angular variation simply 
produces nicely oval contours, with m chosen between 1.3 and 1.5 
for prolate clouds and between 0.7 and 0.8 for oblate ones. We also 
tried using a Gaussian distribution of the form 



71 H , 



c n exp 



+ 



(A4) 



where r and z are as described above for the power law. The pa- 
rameters d n and e n enable the size of the Gaussian (and hence the 
cloud) to be specified. 

The temperature distribution has been assumed to be spheri- 
cally symmetric and varies according to 



T K = citR 03 +b T R~ ' S 



(A5) 



where R is the distance to the centre of the cloud measured in par- 
secs. This distribution, which crudely represents a cloud that may 
be both internally and ex ternally heated, is taken from the results of 
Scoville & Kwan ( 1976) which assumes that (3 — 2. GL98 applied 



this formula to HH24MMS and HH25MMS in LBS23 for mod- 
elling Cl and C O. In principle the temperature variation with 
distance from the centre should be different along the two a xes of 
the flattened model core. However, Scoville & Kwan (1976) show 



that the optical depth probably has a greater effect on the tempera- 
ture fall off than the density law, and thus for simplicity we retain a 
spherically-symmetric temperature distribution. 

For all models a cloud consisting of 60 cylinders and 60 disks 
(i.e. a total of 3540 rings) was used with both the disks and cylin- 



ders packed closer together near the centre of the cloud (this en- 
ables the rapid changes of the power laws near the centre of the 
cloud to be modelled without large jumps in the parameter values 
from one ring to another). The cylinders have an outer radius of 
0.05 pc and the disks extend out to a distance of 0.05 pc above 
and below the centre of the cloud. The exact choice of this dis- 
tance does not affect the results for the Gaussian models. However, 
power law models are affected to varying degrees depending upon 
the steepness of the power law used. For steep power laws the cross- 
section shape remains approximately the same with only a 10 per 
cent change in c n when the radius is increased by a factor of 4. 
For the shallow power law the changes are much bigger. For the 
worst case increasing the radius of the cloud by a factor of 4 leads 
to significantly higher emission at the edge of the cross sections, 
particularly for the 450 /im crossection. This can be corrected by 
making the power law a little steeper (0.1 extra in each direction) 
with a 60 per cent increase in c n . 

In order to model the beam correctly for the rapid changes in 
the centre of the cloud a total of up to 55x55=3025 lines of sight 
are distributed across the beam. This number was determined by 
trying progressively larger numbers of lines of sight for each model 
until the cross-sections stopped changing. 

The program initially calculates the intensity only for the cen- 
tral point in the cloud and then, by adjus ting the values of br (in 
equation \5 ) and c n (in equations 43 & A4 ), it fits the peak tem- 



peratures in the 450 and 850 /im maps. Once the peaks are fitted 
correctly a map is produced and then the four cuts are compared 
with the telescope data. 

With power laws there is usually a problem with deciding 
where to place an inner cut off 7?i n (necessary to prevent a sin- 
gularity at the centre of the cloud). We have used 0.001 pc here, 
although there are only minor changes in the results if this is in- 
creased or decreased by a factor of 5. This is because the extremely 
small area in the centre of the cloud has little affect on the total 
emission over the beam. 

There is also some quest ion as to the precise form of the emis- 
sivity. The Hildebrand ( 1983 ) formula, which we have used here, 
proposes that in the long wavelength regime (A > 250/im) the ab- 
sorption coefficient at a frequency v (GHz) is related to the molec- 
ular hydrogen density iih 2 (cm -3 ) by 



2n Ha 



(— V 

1.2 x 10 25 V750/ 



(A6) 



Others (e.g. Visser et al. 1998 and Testi & Sargent 1998) have re- 
cently suggested that lower values for f3 may be more appropriate. 
We tested this by using the relationship suggested by Testi & Sar- 
gent <|L99§ of «„ = 



9.6x10^ (25o) ' Cm in P laCe ° f e 1 Ua - 

tion|A6l where as before v is in GHz and is in cm -3 . However, 



as mentioned in the text, the precise choice of emissivity law does 
not significantly affect the results (other than the magnitude of the 
temperature and density). 

This paper has been produced using the Royal Astronomical Soci- 
ety/Blackwell Science ETjgX style file. 
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